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ABSTRACT 



_ Non-Fermi liquid behavior is found for the first time in a two-dimensional (2D) system 
Q> ! with non-singular interactions using Haldane's bosonization scheme. The bosonized sys- 



tem is solved exactly by a generalized Bogoliubov transformation. The fermion momentum 



B 

n^ • distribution, calculated using a generalized Mattis-Lieb technique, exhibits a non-universal 

S \ power law in the vicinity of the Fermi surface for intermediate interaction strengths. 
> 



In recent years, there has been great interest in possible non-Fermi-hquid-hke behavior 
in the ground state of two-dimensional (2D) strongly correlated fermionic systems [1]. It is 
well-known that in one dimension the ground state of a fermionic system with arbitrarily 
weak interactions is not a Fermi liquid but a state usually known as a Luttinger liquid [2] . 
For the 2D case, many studies have suggested that, in the limit of weak interactions, this 
behavior does not appear [3]. Studies of systems with a small number of particles also 
suggest Fermi liquid behavior in 2D [4]. Although Khveshchenko et al have found non- 
Fermi-liquid behavior in 2D, they worked only with a singular long-range current-current 
interaction [5]. 

The one-dimensional Luttinger liquid solution can be obtained by a bosonization 
procedure [2]. Haldane has recently generalized this procedure to solve 2D models [6] 
and his method has been further developed in [5]. In this paper, we show using this 
procedure that, even for non-singular short-range interactions, 2D systems may show 
non-Fermi-liquid behavior if the interaction is strong enough. However, in the limit of 
weak interactions it is difficult to distinguish between a Fermi liquid and a non-Fermi 
liquid. 

We diagonalize the bosonized model exactly using a generalized Bogoliubov transfor- 
mation. The fermion momentum distribution n(p) at zero temperature is then calculated 
using a generalization to higher dimensions of the method used by Mattis and Lieb [7] in 
ID. The power law obtained in ID for n{p) near pp, n{p) ~ (p — PfY , also emerges in 
2D. The exponents, 5, can be calculated numerically. We look in particular at the case of 
a 5-function interaction and find that the variation of 5 with coupling constant is similar 
to that found in ID, although 5 is always smaller in 2D than ID for comparable values of 
the coupling constants. 



We begin with a brief description of Haldane's procedure. Consider the fermionic 
system with the spectrum e(p) = f(|p2;| + \Py\)- The spectrum has 4 different branches 
with one in each quadrant of the p-plane. One constructs 4 different density operators 
-Sa(ct) = Y^Ca^p+qCap whcrc a IS the quadrant index and p and p + q he in the a-th 
quadrant. It is easy to show that those operators obey the commutation relations 

[Ba{q), 5/3(-q')] = -<^a/3<5q,q'77r-T^n„ ■ q, (1) 

where A and n^ are the length and normal vector of the Fermi surface in the a-th quadrant. 
Q denotes the volume of the system. In terms of the commutation relations between the 
kinetic energy and Ba{ci), the total bosonized Hamiltonian reads 

^ = ^E(*». + ^)B.(q)B.(-q). (2) 

if the fermionic interaction has the form ^ Z]p,p',qr(p, p', q)Cp_|_qcL Cp/Cp. Since the 
definition of the -Ba(q) excludes contribution to the density operators in which p and 
p + q are in different quadrants, certain scattering terms are not included in the bosonized 
Hamiltonian. These describe the so-called around-the-corner (ATC) scattering. In the 
limit g -^ the contribution of ATC scattering vanishes as q/A. We will return to this 
question later. 

Haldane has pointed out that the 2D system described by (2) will converge to the 
Fermi liquid fixed point in the limit A ^ for non-singular interactions. Haldane's 
observation is obviously correct when the summation over a, (3 is over four quadrants. 
However, if the summation over a and (3 is taken over an infinite number of sectors 
instead of four quadrants, as A ^ 0, the situation may change as we shall see below. 

We consider an interacting fermion model described by the Hamiltonian 

P.« " q p,p',s,s' 



where Planck's constant %, the speed of hght c and lattice length a have been set to unity. 
We assume the spectrum e(p) = f (|p| —pp) where pp is the Fermi momentum. Thus the 
Fermi surface is a circle with radius pp- Following Khveshchenko et al [5], we divide this 
circular Fermi surface into an M-sided polygon (M -^ oo), whose side length is A with 
A <^ j9^. The whole momentum space is thus divided into M sectors. The a-th sector, Tq,, 
is defined by the direction of the unit vector iIq,. In order to understand the low energy 
behavior of the system, we only consider a shell of width L (L ~ A) around the Fermi 
surface. The bosonic operators are defined as -BQ,s(q) = Z]peA„,p+qer„ Ca,p+q,sCaps, where 
a runs from to M — 1. Aq, is the area L x A and centered at p^ , the Fermi momentum 
in the a-th direction. The bosonic operators satisfy the commutation relations of Eq.(l). 

The kinetic energy can be written in terms of the -BQ,s(q) [5] 

(2Tifv 



''kin 



^^^^ q,a,s 



Y,BU^)BU-^) + 0{q^). (4a) 



The interaction is bosonized to give 

+ ^ E V{q)BUci)B,A-^)- 

Here V'{q) = V{q) — Vi{q) with V{q) corresponding to the contribution from direct 
scattering and Vi(g) to that from exchange scattering. For a 5- function interaction, V'{q) 
vanishes, as required by Pauli's principle [8]. In Eq.(3a), we have restricted the momenta 
p and p + q included in the interaction term of (3) to A^ and T^, respectively, as done 
in refs.[6,5]. 

We now solve the bosonized model described hj H = H^in + Hy- The commutation 
relations for the Bas{f\) allow creation and annihilation operators to be defined as follows. 



For iIq, ■ q > 0, 



Basiq) = W77r^^yj4s(-q^ 



(27r 



2 9.' 

(5) 



where Uq^a = I'n.a ■ q|- For n^ ■ q < 0, aj^^(— q) and aas(— q) are replaced by aas(q) and 
^as(q) respectively. The operators aQ,s(±q) obey the canonical commutation relations 
[aas(±q), ajjg, (±q')] = ^c/j^q^q'^ss'. Again we divide the q-plane into M regions. Each of 
them is the region between iIq, and Ua+i, and is denoted by Da- Without losing generality, 
we consider M = AN for convenience. Using Uq^a = ^q,a+2N = ^ci,a+4N, the bosonized 
Hamiltonian can now be rewritten as 

2N-1 

^ = 2E E E<s,q^^,^',qA,s',q. (6) 

n=0 qG-Dn s,s' 

For computational convenience, we organize the 4iV-component vectors, Al^g, so that the 
first N and last N components are the creation operators. For example, we would have 
for the cases n = and n = 1: 

"^0,5 ~ [%sj •■•5 '^Af-l,s5 (^Ns, •••5 03Af-l,S5 %Ar,s5 •••, '^4Af-l,sJ) 

(7) 

"^1,5 — [%,s> •••) '^Af-l,s? (^3N,s, O'N+lyS, ■■■, 0'3n-l,s, O'Ns^ %Af+l,s5 ••■5 '^4Af~l,sJ) 

and so on. Here a^^s — '^ls(~qn) and aa,s = cia,s{q.n) and q„ denotes q G -D„. The 
elements of the matrix Qs,s' are given by 

^«/3,ss'(q) = [(^ + AV'{q))6ss' + AV'(g)5s,_s']c^q,a, for a = P, 
^afs,ss'{q) = AV{q)uJq,lujll]j, for a 7^ /5. 

We introduce the charge density operators, Ra,q, and spin density operators, Sa^q, via 
the relation aQ,(±q) = -^^(-Roq + s Saq)- The Hamiltonian (6) can therefore be expressed 



as H = Hfj + iJp, 

n,D„,a /Q\ 

n,D„,a,/3 

where ^aa„q = {v + AV'{q) ± AV {q))ua,ci with p corresponding to '+' and a to ' — ', and 
^a/3,q ~ '^^ai3,ss' for a ^ {3 . The vectors V correspond to vectors A in (7) with a's replaecd 
by R's. 

The spin density part of the Hamihonian, H„, is aheady diagonal. There are AN 
excitation branches with energies ^aQ,,q • The branches S'^^q and 5'a+2Af,q are degenerate. 

We need to diagonahze the density part of the Hamiltonian, Hp. Formally, this can 
be done by setting up the equation 

/(A) = Det|l]^^+,,,^+,. - ^a\5a^,5^^ = 0, (9) 

where i,j = 1, ..,N, and a,b = 0, ..., 3. With our ordering of the 4A^-coniponents of the 
vector A, Aa = +1 for a = 0, 3 and Aa = — 1 for a = 1, 2. Since Ua^q = uJa+2N,q, we have 
/(A) = /(—A). The 4A^-component 'extended' eigenvectors, X^^\ satisfy: 

QPX(^'^ = A^AX(^\ (10) 

for each A^. Here the matrix A = (AaSabSij). For real A^ [9], we have the extended 
normalization and orthogonality relations 

X-'^'^AX-^ = A', (11) 

where X~^ = {X^^\ ...,X^^'^~^^) is the inverse of the transformation matrix X and A' = 
(A^(5^^/) with A^ equal +1 for A^ > and —1 for A^ < 0. The transformation matrix, 
X = A'X^^ A, can be easily obtained from the orthogonality relations (11). Using 



Eqs.(lO) and (11), Hp is diagonalized to H^ = I] |A^(q)|6'1jq^6'^^q^, where the Oj^^qn ^^^ 
^/i,q„ are true bosonic creation and annihilation operators. 

To determine the nature of the ground state of the Hamihonian (8) and, in particular, 
whether it is a Fermi liquid or not we look at the momentum distribution, n(p). We first 
extend the method used by Mattis and Lieb in ID [7] to higher dimensions. 

A fermion ipsi'^) = Z^pC ^^'^Cps is written in terms of operators with momenta re- 
stricted to each sector, A^,, of the polygonised Fermi surface, i.e. ips{'^) ~ J2a'^as{'^) 
where ipasip) = J2peAa^^^^^^a,p,s- Then the fermion momentum distribution at zero 
temperature is 

n(p) = Y: /rf^P-" < G|C(x)V^.s(0)|G >, (12) 

where jG > is the vacuum state of the interacting system described by the bosonized 
Hamiltonian (8). 

The ground state, \G >, can be related to the non-interacting vacuum state |0 > by 
a canonical transformation of the type \G >= e~*^|0 >, as pointed out by Mattis and 
Lieb. Under this transformation the fermion operator becomes V's(x) = e~*^-?/'s(x)e*^. We 
assume S has the form i J2n,D„,a,i3 ^aQ„^ai3,q„'^f3c[r,- The matrix A can then be determined 
from e^ = X with OajS = AaAa/s- As in the ID case, the fermion operator in the 



bosonic representaion is of the form ■^/'^.^(x) = ^^.^(x)^^"''*^''^ where F{x) 



~ e^p can 



be an arbitrary function of x. J(x) is the current operator and is very analogous to 
its counterpart in ID. Dividing J into J'^ and J^ relating to spin and charge currents 
respectively and solving the 'equation of motion' for f^r = e*'^^-?/'Q,(x)e~*^^ with time 'r' 
[7,5], one obtains 

V^«,(x) = w^,(x)y;(x)^„,. (13) 

Since the elements of the transformation matrix X satisfy Xq,^ = X2N+a,2N+i3, the oper- 



ators Wa and Ya, in the cases ioi a = i = 0, ...,N — 1, are given by 

2-K ,.,-^ Xii — 1 



where D^ = {Dq, ...,Dj^^i+i} and D_ = {D^+i, ...,D2n~i}- Similar relations follow for 
a = N, ...,4N — 1. These results, combined with the orthogonality relations (11), lead, 
after some algebra similar to that in ref.[5], to 

< G|C(x)^as(0)|G >= exp{-Q,(x)} < 0|VL(x)V'as(0)|0 >, (15) 



where 



Qa(x) = — ^ Y. /a(q)(cosq-x-l), 

Y: (X^,(q) + X^ 3^+,.(q)) - 1, if A, > 0, (16) 

j=0 

N-l 

j=0 



/a(q) 



It is seen from (12) and (15) that the momentum distribution n(p) is given in terms of 
the quantities (5a(x). Once these are known we can establish the nature of any singularity 
in n(p) near pp. We note that, although the unobservable quantity /^(q) depends on the 
ordering of {\a} via (16), the physical quantity n(p) does not. 

To illustrate the method, we have solved for the Xq,^ for the case of a simple (5-function 
interaction. In this case V{q) = U and V'{q) = and we can define the dimensionless 
coupling constant g = ttpfU/v. We denote the component of q perpendicular (parallel) to 
Uao for any ag as q± {q\\) and scale them by ty = q\\/L and tj_ = q±/A. For any given M, g 
and t, we denote by fmaxi'^) and Qmaxi'^) the values of /^(t) and (5q(x) corresponding to 
the maximum and minimum values of Aq,, ±Xmax- We also order the Aq, so that A^ = Amax 
for a = ao and A^ = —Xmax for a = ao + 2N. 

8 



We start with some empirical observations on the results we obtain. For small g 
{g < 0.05) we find that fmaxi't) depends on t and is comparable to the other /„. However, 
above g ~ 0.05 there is a crossover to a regime in which /max(t) does not depend on t and is 
significantly larger than the other f^- fa/fmax ~ 10^^ for g = 0.05 while fa/ f max ~ 10~^ 
for g > 1.0. As fmax turns out to correspond to an exponent we denote it by 6. We find 
that, for any given g, 6 converges to a constant as M increases. This is shown in Fig. 1. 
When 6 is much larger than the other fa, we retain only Qmax and set all other Qa to 
zero. This allows us to write 

(17) 



+ Y. j rfMp-k)^(lkl-p^), 



a^ao,ao+2N ' 

where the integration domains, Va, are the sectors corresponding to Da in the full k-space 
and < 0|V'^(x)-?/'(0)|0 >= l/fiX]k^(|k| —pp), the free fermion distribution, is used. 

We now evaluate n(p) from (17). Since the angle a^ can be chosen arbitrarily, it is 
convenient to assign ao in such a way that it labels the sector containing p. The second 
term in (17) then vanishes. Since p G Va^, Qmaxi^) ~ Qmax{x^\,0) and the integration 
over t leaves 

Qmax{xil,0) = 6 ^(l-COSS||), (18) 

Jo S|| 

where W = Lx\\/2. For W :^ 1, this gives Q{x\\,0) = 6 In W"^ (see [5]) and, after 
integrating over x\\ in (18), we obtain 

n(p) = const. — c( )'^sign(|p| —pr), (19) 

ij 

with c = 2^1{5)/7r6 and I{5) = J^ dsn-sT^^ cos S||. The power exponent fmax = S is the 
function of g, shown in Fig. 2. 

For the case of intermediate coupling g > 0.05, where we can obtain the approximate 



form for n(p) given in (19), one can see clearly that the 2D system is not a Fermi liquid. 

In Fig. 2 we compare the variation of 6 with g to the function Sigi) = = — 1 for 

V 1 + 2^1 

the ID case {gi = 2U/v). For the 2D case the non- Fermi-liquid exponent is always smaller 
than in the ID case for comparable values of the dimensionless coupling constants. 

For weak coupling g < 0.05, the contributions from ±Xmax are not dominant and 
contributions from all of Aq, are comparable. We have noticed that the largest of the fa 
varies as lO"^*^""''^) for g ~ 10^". This clearly implies that it will be difficult to distinguish 
between Fermi-liquid and non-Fermi-liquid-like behaviors. This might be the reason why 
the non-Fermi-liquid state has not been found for weak interactions using perturbative 
approaches [3]. 

We would like to comment on the neglect of the around-the-corner(ATC) scattering. 
One might expect that in the limit of a large number of sectors (A ^ 0), ATC scattering 
should become increasingly important, as more and more states are close to the corner of 
a sector. However, as our results show that the exponent, 6, is already well-converged for 
a small number of sectors, we do not believe that neglect of ATC scattering poses serious 
problems. 

Although we have located at the bosonization of the Hamiltonian for a homogeneous 
two-dimensional system, we end up with a Hamiltonian which resembles that for coupled 
chains discribed by Wen [10]. Our results differ from those he obtained for coupled chains 
in one very significant feature. In Wen's case the interchain interactions modify the single 
chain Luttinger liquid state. In our case the interactions between sectors generate the 
non-Fermi liquid state. If we switch off this interaction the system reverts back to a 
Fermi liquid state and not a series of uncoupled ID Luttinger liquids. 

In conclusion, we have discussed a generalization to 2D of the Tomonaga-Luttinger 

10 



model using Haldane's bosonization procedure. The bosonized system has been solved ex- 
actly by a generalized Bogoliubov transformation. The fermion momentum distribution at 
zero temperature has been calculated using an extension of the Mattis-Lieb technique. We 
have found that the momentum distribution exhibits a non-universal power law behavior 
in the vicinity of the Fermi surface for intermediate interaction strengths. It seems to us 
that this is the first time that Luttinger liquid behavior has been found nonperturbatively 
for a two-dimensional model with a non-singular interaction. 
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CAPTIONS OF FIGURES 

Figure 1. The exponents (5 as a function of the number of sectors of the polygonised 
Fermi surface, M, for various values of the dimensionless couphng constant g. The 
exponent, 5, characterises the singularity at the Fermi surface in the momentum 
distribution function, n(p). As M increases, 5 converges to a constant value for 
each coupling constant g. 

Figure 2. The M -^ oo limit for the exponent, (5, as a function of coupling constant g. 
(M is the number of sectors of the polygonized Fermi surface). We also show the 
corresponding results for the ID case. 
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